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ABSTRACT 

We present continent-scale very-long-baseline interferometry (VLBI) - obtained with the 
European VLBI Network (EVN) at a wavelength of 18 cm - of six distant, luminous 
submillimetre-selected galaxies (SMGs). Our images have a synthesized beam width of 
«30milliarcsec FWHM - three orders of magnitude smaller in area than the highest resolu- 
tion Very Large Array (VLA) imaging at this frequency - and are capable of separating radio 
emission from ultra-compact radio cores (associated with active super-massive black holes - 
SMBHs) from that due to starburst activity. Despite targeting compact sources - as judged by 
earlier observations with the VLA and MERLIN - we identify ultra-compact cores in only 
two of our targets. This suggests that the radio emission from SMGs is produced primarily 
on larger scales than those probed by the EVN, and therefore is generated by star formation 
rather than an AGN - a result consistent with other methods used to identify the presence of 
SMBHs in these systems. 

Key words: galaxies: starburst - galaxies: formation - cosmology: observations - cosmol- 
ogy: early Universe 



1 INTRODUCTION 

With bolometric luminosities rivalling quasars, submillimetre 
galaxies (SMGs) are some of the most extr eme objects in the 
Universe (for a review, see iBlain et al.l [2002t> . yet they outnum- 
ber similarly luminou s quasars by many orders of magnitude 
( Chap man et alj [2005h . Selected via their dust-reprocessed ther- 
mal emission, they are thought to be powered by intense bursts 
of star formation triggered by major mergers of gas- rich galaxies 
Tacconi et alj200dJ2008l ; lYounger et alj2008all2oTch at z ~ 2-3 
Chapman et alj2005l) with a significant tail extending out to higher 



redshifts teales et al.l200ilYounger et alj200ll2009al:IWang et al 



2007 



2008 



20091: iGreve et al.ll2008l: ICapak et alJl2008l : ISchinnerer et al 
Daddi et al. 2009; ICoppin et al.l2009n . and are the likely pro- 



genitors of massive galaxies in the local Universe iScott et alj2002l : 
IBlain et alj|20ol : ISmail et alj|2004l : ISwinbank et alj|2006h . In ad- 
dition, it is thought that infrared- (IR-)luminous objects (Lir Z 
10 12 ~ 13 Lq) come to dominate th e cosmic star-formation rate den- 
sity a t high redshift (z ~ 2 - 4; B lain et al.ll 19991 ; iHopkins etal] 
l201Ct Hopkins & Hernquistl l2010h . making the SMG population 
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an important contributor to the build-up of stellar mass during the 
epoch of galaxy formation. 

At the same time, over the past twenty years it has become 
clear that super-massive black holes (SMBHs) in the nuclei of 
galaxies are common, even at high redshift (e.g. llvison et alj2008h 
and that their masses are strongly correlated with the stellar com- 
ponent of their host g a laxies (e.g., iKormendv & Richstonelll99: 



Mag orrianetal.il 1998 
2000; iTremaine et al 



Gebhardt et al. 2000; Ferrarese & Merritt 
20021 ; iNovak et alj I2OO6I : IHopkins etal 



2007b). These correlations are indicative of a close and fundamen- 
tal link between SMBHs and stellar populations, and the enormous 
difference (~ 10 9 ) in linear scales suggests that this may be accom- 
plished via accretion-related radio jets or radiatively-driven winds. 
Recent theoretical models cast the co-evolution of SMBHs and 
galaxies in the contex t of a cosmic cycle driven by major mergers 
( IHopkins et alj|2008bl lal and references within), wherein gas-rich, 
spiral galaxies collid e and trigger a %2& inflow, thereby fuelling 
a nuc lear starburst jHernquistT 1989; Mihos & Hernquist 1 19941 
199fJ) and feedba c k-regulated SMBH growth fsllk & Reesll 199^ : 
Page et alj |2004 |Pi Matteo et alj |20O5t IHopkins et alj l2007iJ : 



Younger et a l. 2008b), eventually revealing a bright quasar, after 



which gas exhaustion and v iolent relaxation transform the remnant 
into a red elliptical galaxy (Barnes &Hernquistlll992h . 

These models of galaxy evolution, which match constraints 
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from both galaxy Jllopkins et alj 20061) and SMBH popula- 
tions |Robertson et alj|2006l : iHopkins et alj|2007al ; lounger et al.1 
l2008bn T predict that hyperluminous starbursts at high redshift will 
be associated with periods of active SMBH growth. SMGs are 
thus likely candidates for the transition objects that theory pre- 
dicts should be powered by a m ixture of star formatio n and ac- 
tive black hole growth (see also lYounger et al. 2009b). To date, 
searches for actively gr owing SMBHs in SMGs ha ve focused 
on X-ray observations dAlexander et al] l200j |2008h and mid- 
IR photometry and spectroscopy Uvison et alj 120041: iLutz et alj 
20051; iMenendez-Delmestre et al.ll200ll2009l ; IValiante et al.ll2007F 
Pope etalj 12008). Whilst promising, these approaches are often 



Table 1. Names of SMGs targeted by our EVN observations. The six are 
split into two subsamples, the first three being significantly brighter. 



compromised: even hard X-ray searches will miss Compton-thick 
SMBHs, and the interpretation of mid-IR spectra can be very sen- 
sitive to modelling assumptions, especially regarding the origin of 
power-law spectra dYun et afl200ll ; lYounger et alj|2009bl) . 

Very-high-resolution radio imaging, provided by continent- 
scale very-long-baseline interferometry (VLBI), avoids many of 
these issues; even in the most extreme environments, the interstellar 
medium is optically thin to radio emission and the dominant radio 
emission mechanisms are well understood. In fact, there is a well- 
known upper limit to the b rightness temperatu re (Tb) a starburst 
can achieve (Tb < 10 5 K - ICondon et alj|l99ll) , and any compact 
object detected at VLBI resolution will exceed (or be close to) this 
limit. 

A number of authors have used this technique in order to dis- 
tinguis h between AGN and starburst activity in both high-redshift 
QSOs ( iBeelen et al.l2004l:lMomjian et alfeOOi koOV) and in ultra- 
luminous infrared galaxies (ULIRGs), the low-redshift analogues 
of the SMGs. The latter are particularly important as they present 
a far more detailed view of the relationship between the star for- 
mation and AGN emission in an intense starburst than will ever be 
possible in a distant SMG, and aid in the interpretation of these 
less spatially resolved sources; the synthesised beam of a typical 
VLBI array is high enough, for example, to resolv e individual su- 
perno va remnants in local U LIRGs e.g. Arp 220 ([Lonsdale et all 
l200fj) and IRAS 17208-0014 TMomiian et alj J2006T) . Examples of 
ULIR Gs that are believ ed to contain an AGN include Arp 220 
dDownes & Eckartll2007l) and Mrk 273 JCarilli & TayloihoOOl) . 

Only one VLB I observation of a c lassical, high-redshift, 
optically-faint SMG (Momi ian et al ] |20ld) has been published to 
date, continuum emission from GOODS 850-3 having been de- 
tected in a tapered High Sensitivity Array (HSA) image. In this 
paper we present the first VLBI survey of SMGs, using the Euro- 
pean VLBI Network (EVN) to search for ultra-compact cores in a 
sample of objects in the Lockman Hole. These data, with a synthe- 
sized beam size of ~30milliarcsec (mas) FWHM, are the highest- 
resolution radio detections of SMGs ever achieved and have al- 
lowed us to put the tightest constraints to date on the brightness 
temperatures of a significant sample of high-redshift starbursts. 

This paper is organised as follows: in |j2]we describe our target 
selection and in ijUgive details of the observing strategy, data re- 
duction and imaging. In Sj4]we present our high-resolution images 
of each SMG and describe each source in detail. In ij5]we discuss 
our findings before presenting our conclusions in Sj6] 

Where necessary we have assumed a flat ACDM cosmol- 
ogy of n A = 0.73, n m = 0.27 and H = 71 kms^Mpc -1 
(Hinsh awetalj2009l) . 
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2 TARGET SELECTION AND AUXILIARY DATA 



The six targets were drawn from the Biggs & Ivison (2008) sam- 
ple of SMGs in the Lockman Hole with high-resolution (synthe- 
sized beam size ~200-500mas) radio imaging provided by the 
Multi-Element Radio-Linked Interferometer Network (MERLIN), 
an array of seven telescopes across the United Kingdom. Table Q] 
gives the various names of each SMG (most are detected in sev- 
eral submm surveys at different wavelengths) whilst in Table [2] we 
reproduce various observational data, including the position of the 
radio counterpart, their sub mm and 1.4-GHz flu x densities (mea- 
sured both by MERLIN -iBiggs & Ivisor]|2008l - and the VLA - 
llvison et"aT]|2002 ; IBiggs & IvisorJ 2006h and radio spectral index 
l Ibar et al. 2009L 20 id) . A spectroscopic redshift is given where 
available. The first three SMGs were selected on the basis of their 
high peak flux density in the MERLIN map, as well as the presence 
of compact emission that could not be resolved at full MERLIN 
resolution. In addition, with a maximum separation of 5.7 arcmin, 
they could be observed with a single pointing of the EVN (the pri- 
mary beam of the Effelsberg telescope is 6 arcmin at 18 cm; see 
Fig-UJ. In order to widen the range of luminosities in the sample - 
and hence to probe galaxies representative of the wider SMG pop- 
ulation - the second group of SMGs are considerably fainter (by 
3 x , based on their VLA flux densities) and also detected by MER- 
LIN. The maximum separation between these three sources is only 
2.2 arcmin, so the primary beam attenuation is reduced compared 
to EB038 - useful when trying to detect such faint sources. 

The 850-^im data for this sample was provided by the Sub- 
millimeter Common Us er Bolometer (SCUBA - iHolland et alj 
ll999l : ICoppin et al]2006|) on the 15-m James Clerk Maxwell Tele- 
scope, the 350-^im data are from SHARC-II imaging at the 10- 
m Caltech Submill imeter Observatory (CSO - iKovacs etalfeoOtjl; 
ICoppin et al]2008l), and the millimetre (A 1100 - 1 200/xm) sur- 
vey data dGreve et al 1200 4l; lLaurent et alj |2005L|2006|) are from the 
MAMBO bolometer dKrevsa et alj|l998h on the Institut de Radio 
Astronom ie Millimetrique 's (IRAM's) 30-m telescope and from 
Bolocam faaig et alj|2004l) on the CSO. 



3 OBSERVATIONS, DATA REDUCTION AND IMAGING 

Our EVN observations of the Lockman Hole were split into two 
separate observing sessions, the first targeting the brightest three 
SMGs (project code EB038; 2008 June 3 and 5) and the second 
(EB041; 2009 June 5 and 6) targeting the remaining three (Fig. [TJ. 
The EB038 and EB041 observations had durations of 12 and 14 hr. 
The seven telescopes which participated in EB038 are detailed in 
Tableland include the 76-m Lovell telescope at Jodrell Bank, the 
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Figure 1. VLA 1.4-GHz map of the Lockman Hole showing the locations of the six targeted SMGs and their division into two separate observing epochs. The 
dashed circles correspond to the primary beam of the Effelsberg telescope at 1 8 cm centred on the pointing positions of all telescopes apart from the Westerbork 
Synthesis Radio Telescope (WSRT) which pointed at each SMG consecutively; the small green boxes around each SMG position are approximately the size 
of the WSRT field of view (11 arcsec). The pointing positions (red crosses) for all telescopes (apart from WSRT) were 10 h 52 m 01! 250, 57° 18' 42'.'76 for 
EB038 and 10 h 51 m 58?150, 57° 26' 01'.' 13 forEB041. The triangle just outside theEB041 field marks the source that was detected using the SMG02 data. 



100-m telescope at Effelsberg and the Westerbork array of 14 x 25- 
m antennas. EB041 also included the 25-m Urumqi dish in Western 
China, this providing much longer baselines and therefore tighter 
constraints on the brightness temperature, 7b. Due to a mechanical 
failure of the Lovell telescope shortly after the beginning of EB041, 
this was replaced by the 25-m Mark II telescope (also located at 
Jodrell Bank) for the vast majority of the project. 

The faintness of the targets made phase referencing essential 
and for our calibrator we used the unresolved ~200-mJy source, 
J1058+5628, which is located only one degree from the SMGs. For 
EB038 we used a cycling time of 300 s, but increased this to 600 s 



for EB041. As explained in S|2] the targets lie within the primary 
beams of the largest telescope and so could be observed in a single 
pointing. An exception was the WSRT which in phased-array mode 
has a field of view of only 1 1 arcsec. Therefore, this telescope was 
cycled consecutively through the three SMG positions, changing 
after every second SMG scan. This strategy required that the data 
be correlated three times, once at each of the three WSRT pointings. 
The bright quasars 3C 345 and 4C 39.25 were both observed for 
600 s as bandpass and delay calibrators. 

Both EB038 and EB041 were observed with dual circularly- 
polarised receivers tuned to a wavelength close to 18 cm. The to- 
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Table 2. Observa t ional data of our sample. 850-/im data are ta ken fromlCoppin et alj f2006l) , 1100-/im data from lLauren t et al. 12005) and those at 1200 
from lGreve et 350-/xm flux densities are taken from lLaurent et alj j2006f) . The positions were measured from VLA 1 .4-GHz imaging llvison et alj 

l2002l;lBiggs & Ivisonll200d) and all redshifts are spectroscopic. The VLA total flux densities and the MERLIN peak flux densities have both been measured 
at 1.4 GHz and corrected for bandwidth smearing and primary beam attenuation. Due to the different antenna sizes, for the MERLIN data the attenuation 
has been calculated using a sensitivity-weighted average of each baseline. The MERLIN maps can display multiple components: the quoted fluxes refer to 
the brightest. The fina l column gives the spectral index between 1.4GHz (VLA) and 610 MHz (from imaging with the Giant Metre-wave Radio Telescope - 
Ilbaretall200ll2010l) where S v oc u a . 



Name 


Q J2000 
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<5j2000 
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(mJy) 


(mJy) 


Sll00|120UAtm 

(mJy) 
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Si. 4, VLA 


•Si. 4, MERLIN 
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1 4 
"0.61 


SMG06 


10:52:07.49 


+57:19:04.0 


38.0 ±7.2 


4 7+1-5 
' -1.6 




2.689 ± 0.01 


246 ±8 


296 ± 8 


±0.83 


SMG10 


10:51:41.43 


+57:19:51.9 






4.1 ±0.9 


1.212 ± 0.01 


295 ±9 


157 ±10 


-0.69 


SMG11 


10:52:13.38 


+57:16:05.4 






6.0 ± 1.4 




245 ± 5 


209 ± 11 


-0.66 


SMGOl 
SMG02 
SMG04 


10:52:01.25 
10:52:04.23 
10:51:51.69 


+57:24:45.8 
+57:26:55.5 
+57:26:36.1 


24.1 ±5.5 
24.9 ±9.1 
25.7 ± 15.8 


8.8 ± 1.0 
10.6+^ 


3.6 ± 0.6, 4.4 ± 1.3 
3.6 ±0.6 
1.6 ±1.6 


3.38 ± 0.02 
1.482 ± 0.01 
1.147 ±0.01 


73 ± 5 
66 ± 5 
110 ±9 


60 ±9 
64 ± 10 
76 ±11 


±0.21 
-0.56 
-0.43 



Table 3. Details of the telescopes used in the EVN observing sessions in 
2008 (EB038) and 2009 (EB041). The Lovell telescope suffered a mechan- 
ical failure shortly after beginning EB041 and only contributed about an 
hour of data before being replaced by the smaller Mark II telescope. 



Telescope Diameter (m) Project 



Lovell, UK 


76 


EB038 


Mark II, UK 


25 


EB041 


Effelsberg, Germany 


100 


both 


Westerbork, The Netherlands 


14 x 25 


both 


Onsala, Sweden 


25 


both 


Medicina, Italy 


32 


both 


Noto, Italy 


32 


both 


Torun, Poland 


32 


both 


Urumqi, China 


25 


EB041 



tal bandwidth was 64 MHz (8 x 8-MHz bands) sampled at 2 bits 
for a total data rate of 512 Mbits s . The data were correlated at 
the Joint Institute for VLBI in Europe (JIVE) with the EVN Mark 
IV Data Processor; bandwidth and time smearing losses were re- 
stricted to ~10 per cent at a radius of ~1 arcmin from the phase 
centre of each image by averaging the correlator output every 1 s 
and splitting each 8-MHz band into 64 channels. Unfortunately, 
during the correlation of data for SMG02, the disk containing the 
data from the Effelsberg telescope failed. As a result, only one quar- 
ter of its data were salvaged and the sensitivity of the images of this 
source and of SMG04 are significantly worse than those of the other 
targets. 

The data for each target/epoch were reduced using AIT'S in 
the following way: initial amplitude calibration was accomplished 
by applying system temperatures measured during the observa- 
tions, with the standard gain curves of each antenna. After finding 
and removing corrupted data (as well as the WSRT scans corre- 
sponding to the other two SMG positions) the delay error on each 
antenna was calculated from a small portion of 4C 39.25 data and 
removed from the entire dataset. The same source was also used 
to measure and flatten the bandpass shape in both amplitude and 
phase. Residual delay, rate and phase solutions were then calculated 
for the phase calibrator and interpolated onto the target SMG. Fi- 
nally, time-dependent amplitude gain corrections were determined 
for J1058+5628 and applied to the target data. 



Once calibrated, the target data were imaged using IMAGR and 
naturally weighted to ensure maximum sensitivity. For the EB038 
data, the combination of this and the uv coverage resulted in a 
nearly circular synthesised beam of 27 x 25 mas 2 and an r.m.s. 
noise level of 10.5 /iJy beam - 1 . For the EB041 data, we found that 
applying a Gaussian taper to the weights (with a value of 30 per cent 
at a distance of 5 MA in both the u and v coordinates) produced 
the best images. This produces r.m.s. noise levels between 10.7 
and 14.3 /ijybeam -1 and a synthesised beam with dimensions of 
23 x 22 mas 2 . 



4 RESULTS 

The final images are shown in Fig. [2] Of the six SMGs that we tar- 
geted, two (SMG06 and SMG1 1) have been detected. Of the other 
four, we have searched for emission on larger scales by tapering 
the data in the (u, v) plane. The shortest baseline in our EVN array 
(between the WSRT and Effelsberg) is also the most sensitive and 
corresponds to emission on scales of approximately 100 mas, still a 
factor of four smaller in beam area than that of the MERLIN maps. 
This, combined with the tapered maps having sensitivities that are 
poorer by about 50 per cent and the higher sidelobes in the syn- 
thesised beam, resulted in the continued non-detection of all SMGs 
other than SMG06 and 1 1. 

The lack of detections from the EB041 data might be ascribed 
to problems with the observations (especially the phase referenc- 
ing) and/or the data reduction. However, as our data were correlated 
with a large number of channels and short integration times we can 
search for emission from sources elsewhere in each SMG's field. 
Three clear detections were made, one in each field and are shown 
in Fig. [3] We can be confident therefore that the phase referencing 
was successful and that the failure to detect emission from SMGOl, 
SMG02 and SMG04 is due to the lack of suitably-bright compact 
emission in these sources. 

For the two detections, we have measured their flux densities 
using JMFIT in AIT'S. The fitted flux density values need to be 
corrected for the attenuation of the primary beam and, as with the 
MERLIN flux densities, this is not straightforward because of the 
different diameters of the telescopes. We have therefore taken the 
weighted mean of the attenuation factors for each baseline, assum- 
ing that the primary beam of each is Gaussian in form and using the 
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Figure 2. Images of our sample of SMGs. LEFT COLUMN: robustly-weighted VLA images at 21 cm, CENTRE COLUMN: naturally-weighted MERLIN images 
at 21 cm, RIGHT COLUMN: naturally-weighted EVN images at 18 cm. For the VLA maps, contours are plotted at —3, 3, 6, 12, 24, 48x the r.m.s. noise. The 
MERLIN and EVN maps have contours at —3, 3, 4, 5, 10, 20, 30, 40x the r.m.s. noise (Table|4)- Synthesised beams are shown in the bottom-left corners and 
each image is the same multiple of the beam area. Note that the MERLIN and EVN images shown here have not been corrected for primary-beam attenuation. 



IStromld2004h formulism for the primary beam of an interferometer 
pair. The resulting values and the corrected peak flux densities are 
both shown in Table [4] 

The attenuation factors are generally insignificant, apart from 
SMG10 and SMG11. As already noted, all MERLIN flux densi- 
ties have been corrected for bandwidth smearing and primary beam 
attenuation, the latter calculated in a similar fashion to the EVN 
values. The pointing centre of the MERLIN observations was the 
centre of Fig. Q] and so the attenuation factors are in general larger 
than for the EVN. 

We now give brief descriptions of each source: 

SMG01: The first source of the faint subsample is actually the 



brightest submm source in the Lockman Hole East region (850.01). 
Recent Spitzer mi d-IR spectroscopy h as determined a reliable red- 
shift, 3.38 ± 0.02 JCoppin et alj2010l) . SMG01 is detected as a sin- 
gle component by MERLIN with a 1 .4-GHz flux of 90 /iJy, a mor- 
phology that is consistent with its flat spectrum, a = +0.21. This 
would have corresponded to a 9-a detection were it unresolved at 
the EVN resolution, but only weak (3-cr i.e. ~30 /iJy beam -1 peak) 
emission is seen near the centre of the map. 

SMG02: The MERLIN map shows two distinct components, but 
neither are detected in the EVN image. 

SMG04: The MERLIN map is dominated by a single component 
which is elongated slightly, north-south. There is no significant de- 
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Figure 2 - continued 



tection in the VLBI map, although we note that there is weak (3-cr 
i.e. ~45 /ijybeam" 1 peak) emission in the EVN map close to the 
expected position of the source. 

SMG06: This source is detected robustly by the EVN, is found to 
be unresolved by JMFIT and has a flux density that is significantly in 
excess of the VLA value, presum ably due to the source variability 
first noted bv llvisonet all (2002). The MERLIN image is slightly 
resolved, but has a peak flux density (296 /iJy beam -1 ) that is simi- 
lar to the EVN value. SMG0 6 has a highly-invert ed radio spectrum, 
a = +0.83, as also noted bv llvison et aT]b002h . 
SMG10: Despite being the brightest radio source in our sample 
(VLA flux density = 295 /iJy), SMG10 is not detected by our EVN 
observations. An examination of the MERLIN map shows that this 
source contains two components, although the brighter of the two 
dominates, with a peak flux density of 157 /^Jybeam -1 . It is pos- 
sible that the weaker of the two components represents a jet, but it 



could equally be one component of a late-stage merger. The radio 
spectral index of the entire system, a — —0.69, is consistent wit h 
an optically- thin radio jet or a ca nonical starburst jlbar et al.l2o"Toh . 
As noted bv llvison et all (2007), matched high-resolution imaging 
at multiple frequencies are required to disentangle the various emis- 
sion components in sources like SMG10. 

SMG11: This source is also detected by the EVN, but at a consid- 
erably lower significance than SMG06, partly due to the primary- 
beam attenuation. Even taking this into account, only 38 per cent 
of the VLA flux is recovered. The MERLIN image is dominated 
by a single (resolved) component with a peak flux density of 
209/iJybeam -1 . The EVN detection is also marginally resolved, 
approximately along the north-south axis. The deconvolved size 
along this axis is 20 mas which is equivalent to a spatial extent 
of 165 pc at the assumed redshift of 2.4 for this source. 
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Figure 3. EVN 18-cm maps of VLA- and MERLIN-detected sources that lay within the primary beam of the EB041 observations. Two of the sources are 
included in the control sample o f Bigg s & Ivisonl (2008) and are labell ed as such in the figures - C06 was imaged from the SMG01 data and C08 from the 
SMG04 data. The source in the SMG02 field was not included in fhe lBiggs& Ivisorj 12008) sample, but is marked as a triangle near the top of Fig. [T] These 
detections demonstrate that our observational and data reduction procedure is sound i.e. the lack of detections of SMGs from the faint sample is due to their 
faintness. The synthesised beams in these images have dimensions of approximately 24 X 20 mas 2 and contours are plotted at —3, 3, 4, 5, 7, 10, 20x the 
r.m.s. noise (16-20 /Jy beam -1 ). Data from the Urumqi telescope are not included so as to minimise the effects of smearing. The WSRT data are also not 
included as all sources are outside its primary beam. 



5 DISCUSSION 

The power of true, continent-scale VLBI instruments like the EVN 
to identify actively growing SMBHs lies in their ability to dis- 
tinguish between different emission mechanisms via Tb, sepa- 
rating ultra-compact r adio cores from mo re extended emission 
dLonsdale et~ai1ll993l) . ICondon et all Jl99lh showed that there ex- 
ists a maximum brightness temperature, Tb, for radio emission as- 
sociated with a starburst. Assuming a canonical starburst i.e. syn- 
chrotron, with a spectral index of —0.8, an electron kinetic tem- 
perature of 10 4 K and scaling the observed temperature to the rest- 
frame of the SMG, we derive maximum values of Tb of approxi- 
mately 5 x 10 4 K for all our sources. Both of our EVN detections 
shown in Fig.|2]have Tb Z 10 5 K (Table[4}, significantly above this 
critical value and both are therefore extremely likely to be produced 
by an active SMBH. 

For the SMGs that are not detected, Table|4]shows that at most 
one half of the radio emission detected by the VLA could be con- 
tained within a compact - i.e. VLBI - component and therefore that 
the implied upper limits on Tb are rather low, < 10 5 K for each. 
Although the limits are not quite low enough to definitively rule 
out the presence of a radio-emitting AGN in these SMGs, if such a 
component is present it must be very near the level of detectability 
in our EVN images. 

What is clear though is that for these four systems the bulk 
of the radio emission must arise on scales larger than the EVN 
beam. This can be attributed to several mechanisms: cosmic rays 
accelerated by supernova remnants created by short-lived, mas- 
sive star s (i.e. star-form i ng galaxies obeying the far-IR/rad i o cor- 
relation: ICondorj 1 19921: [Thompson et ail l200d Rbar et al.1 120081 : 
lLacki et alj2009l : lLacki & Thompsonl2009l . llvison et al.l2010ah . or 
low-surface-brightness AGN-driven radio emission, e.g. jets. Given 
that the presence of an AGN seems unlikely given the non-detection 
of compact components in 4 out of 6 of the SMGs, the evidence 
favours the first explanation, i.e. that the radio emission in these ob- 
jects is predominantly produced by a starburst. Furthermore, since 
this sample was selected for compact emission at ~300-mas reso- 
lution and as most SMGs display extended radio emission at that 



resolution (Biggs & Ivison 2008), these results suggest that the ra- 
dio luminosity produced by SMGs is generally dominated b y star 
formation, a similar finding to that of iMomiian et al. who, 
in the only other published VLBI observation of an SMG to date, 
find that GOODS 850-3 is also consistent with star formation being 
the origin of the radio emission. 

It is also possible that the clearest case of an AGN detection in 
our sample, SMG06, may not actually be a starburst galaxy at all, 
but a highly-obscured AGN-dominated system. The presence of an 
AGN is virtually guaranteed by the source variability that results in 
the EVN flux being in excess of that detected by the VLA, but the 
large EVN flux also leaves very little or no flux that might originate 
from a powerful starburst. In this scenario, the slight resolution of 
the MERLIN image might be due to a radio jet that is resolved out 
by the EVN beam although we stress that the evidence for resolu- 
tion in the MERLIN image is slight, the total flux only exceeding 
the peak flux by 9 per cent. However, given that the source has 
clearly brightened since the VLA observations in 2001, simultane- 
ous EVN and VLA imaging might have detected a clear excess of 
flux with the larger beam. 

That starburst emission is dominant in our sample is not 
necessarily unexpected, given the prediction s of current mod- 
els of SMG formation and evolution (e.g . . iBaugh et alj 120051: 



I Swinbank et"aT]|2008l : iNarayanan et alj|2010l l2009allbl; 



. iDav e et al 

2010) and observations (e.g.lChapman et al.l2004l:lAlexander et al 
2003 : IValiante et al.|[2007l: I Younger et alj|2008al, bold: iPope et al 
| 2008| : iBiggs & Ivisonl 120081 : iMenendez-Delmestre et al.1 12001 
1 20091 : llvison et al.1 l2010bl) which suggest that the luminosity of 
SMGs is dominated by star formation. However, it is not neces- 
sarily consistent with observations of the far-IR/radio correlation 
in these systems. While several autho rs have claimed that the local 
corre l ation holds at high redshift (e.g.lGarrettl2002l:lAppleton et al.1 
l2004l : llbar et al]2008l : lYounger et al.l2009cl : lMurphvll2009h . SMGs 
have been found to fall systematically below the local relation. 
This is usually taken to be indicative of AGN activity (e.g. lCondonl 
Il992h and therefore potentially in conflict with a diagnosis of star- 
formation-dominated radio emission. 

In Table [5] we summarise the far-IR/radio properties for the 
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Table 4. Measured noise levels, flux densities and brightness temperatures for our six targets. Where the target was not detected, we give a 3-cr upper limit. 
Column 3 contains the primary beam attenuation which is used to correct the values given in the following columns. Brightness temperatures have been 
corrected for redshi ft (T b = (1 + z)T b obs ) using the spectroscopic redshift except for SMG1 1 where we have assumed the median SMG redshift of z = 2.4 
I Chap man et alj2005t) . These redshifts have also been used to calculate I Wright 2006) the spatial extent of the beam at the distance of the source. 



Name 


r.m.s. noise 


Primary beam 


EVN peak flux 


EVN total flux 


Recovered VLA 


T b 


Spatial Resolution 




(/jjy beam ) 


attenuation 


(/xjybeam ) 


GuJy) 


flux fraction 


(K) 


(pc 2 ) 


SMG06 


10.5 


0.98 


306 ± 11 


306 ± 11 


1.24 


5.0 x 10 5 


218 x 197 


SMG10 


10.3 


0.85 


< 37 




< 0.13 


< 3.6 x 10 4 


227 x 205 


SMG 11 


10.6 


0.83 


83 ± 13 


92 ± 23 


0.38 


1.4 x 10 5 


224 x 203 


SMG01 


10.7 


0.97 


< 33 




< 0.45 


< 8.5 x 10 4 


181 x 175 


SMG02 


11.7 


0.97 


< 36 




< 0.55 


< 5.2 x 10 4 


200 x 188 


SMG04 


14.3 


0.98 


< 44 




< 0.40 


< 5.5 x 10 4 


188 x 184 



sample presented in this work. All avail able far-IR photo metry 
was fitted with an isothermal greybody (Hildebrand 1983) with 
Sv ~ B u (Ta) - where /3 is the dust emissivity, Td is the physi- 
cal dust temperature, and B v is the Planck function - and the errors 
in the fitted parameters were estimated using a Monte Carlo simu- 
lation assuming Gaussian errors in the photometric data. This form 
of the greybody function assumes the interstellar medium is op- 
tically thin in the far-I R, which may break down in the most lu- 
minous starbursts (e.g.. [Scoville et al] Il99ll : iDownes & Solomon! 
1 19981 : IPapadopoulos et aL 2010h . butwhile the fitted temperature 
is sensitive to the the full emissivity term, Ltr and q (the loga- 
rithmic far-IR/radio ratio) are robust for v$ > 1000 GHz (see also 
lYun & Carillil2002h . 

For SMG10 and SMG11, which have only one photometric 
point avail able, we have fixed Td = 35 K, app ropriate for typi- 
cal SMGs dKovacs et al.l2006l : ICoppin et al.l2008h . and for SMG1 1 
we have assumed z = 2.2+qI; (the med ian and interquartile range 
for SMGs with spectroscopic redshifts - Chapm an et al.l2005T) . We 
have also fixed j3 = 1.5, which is consis tent with local ULIRGs 
dDunne etalJfe OQO; Yang & Phillipsl l2007u as well as current dust 
models (Weingartner & Draine 2001). For the radio, we adopt ra- 
dio spec tral slopes derived from 610-MHz maps obtained with the 
GMRT dlbar et alj|2009l |2010h in conjunction with the VLA 1.4- 
GHz data. As with the mean q estimated for SM Gs with 350-/im 
photometry (q = 2.09 ± 0.09: iKovacs et alJ|2006T) . the objects in 
this sample fall below the local value. In fact, the mean for the 
star-formation-dominated SMGs with no detection in EVN imag- 
ing is (q) = 1.8 8 ± 0.22; well below the local value of q = 2.34 
dYunetal.ll200lh . 

There are three possible explanations for this discrepancy. 
First, owing to a combination of sparse sampling and bad approx- 
imati ons (e.g. an isothermal dust population; see IClements et all 
l2010h . it may be that the far-IR luminosity is systematically un- 
derestimated and these objects do, in fact, fall on the local rela- 
tion. S econd, despite some observational evidence to the co ntrary 
(e.g.. lGarrettl2002l ; lAppleton et alj2004l ; lYounger et al.l2009ch , the 
far-IR /radio correlation co uld evolve with redshift and/or environ- 
ment dlvison et alj|2010ah . Finally, because SMGs are thought to 
be driven by major mergers ( Greve et al. 2005; Tacconi et al. 2006, 
3LI2O 



120081 : lYounger et al.ll2008al . l2010l : iNaravanan et alJl201()lT 2009a), 
there could be luminous synchrotron bridges analog ous to those 
seen i n interacting systems in the local Universe Condon et al.l 
1 19931 120021 ; iMurphv 2009). We note that there are also several 
more prosaic explanations, e.g. sample bias, flux boosting and/or 
mis-identifications in the radio waveband. 

For those targets with detections - SMG06 and SMG 11 - 
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Figure 4. SMBH mass versus stellar mass for all six SM Gs, with black 
hole masses estimated from the 5-GHz radio luminosity jLacv et al]|200lh 
and stellar masses estimated from fitting SED templat es to the optical-IR 
photometry iDve et alj20O8l) . The local relation from Haring & Rix] d2004l) 
is denoted by the solid line and its intrinsic scatter as d ashed lines. Errors 
on the black hole masses are equal to the scatter in the lLacv et all d200ll) 
relationship (1.1 dex) and those on the stellar masses are calculated as doc- 
umented in lDve et all fe008l) . At least one of the EVN non-detected SMGs 
is not consistent with the local relation and, due to being upper limits, all 
of the non-detections might signify black holes that are undermassive com- 
pared to the nearby universe. 



we can use their inferred radio luminosity to estimate the SMBH 
mass. Radio observations of quasars - including both radio-loud 
and radio-quiet objects - have established a stro ng correlation be - 
tween their radio luminosity and SMBH mass dLacv et alj[200lh . 
Assuming that this relation holds for AGN in heavily obscured en- 
vironments, and that accretion is Eddington-limited (which is ex- 
pected during the peak o f both starburst and SMBH activity from 
merger-driven models - ISpringel et aT] |2005). then this observed 
relation is: 

log(Af BH /M ) = 0.52 x log(L 5 GHz/ W Hz" 1 sr" 1 ) - 4.2 

where Tsghz is measured in the rest frame. Though the intrinsic 
scatter is large - approximately 1.1 dex in Mbh at fixed Tsghz - 
we can use these observations to provide a rough estimate of the 
masses of SMBHs detected in the EVN imaging. Again, for the 
radio we adopt the spectr al slopes determined us ing 610-MHz maps 
obtained with the GMRT Jlbar et alj2009ll2010t) and estimate black 
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Table 5. Derived dust temperatures (T d ), luminosities and far-IR/radio ratio (q) for each SMG. We define the far-IR luminosity between rest-frame 40 and 
1000 /im. For SMGs 10 and 11, only one photometric point was a vailable in the far-IR and so we fixed to 35 K, consistent with the mean SMG far-IR 
spectral energy distribution inferred from 350-1200 fim photometry iKovacs et al. 2006; Coppin et al. 2008). As a spectrosc opic redshift was not available for 
SMG 1 1, we have assumed z = 2.2^'g (the median and interquartile range for radio-detected SMGs- Chapman et alj 2005) and propagate those uncertainties 
into the quoted errors on the fitted parameters. Also shown are the derived stellar masses, calculated following Dye et al. 1 2008) and plotted in Fig. [4] 



Name 


T d 


log(L m /L ) 


q 


log(M m /M Q ) 




(K) 








SMG06 


56 ± 15 


13.1 ±0.2 


2.32 ±0.25 


11.6 ±0.2 


SMG 10 


35 


12.7 ±0.1 


1.93 ± 0.12 


11.7±0.2 


SMG 11 


35 


12.8 ± 0.1 


1.45 ± 0.35 


12.7 ±0.1 


SMG01 


28 ±2 


12.8 ±0.1 


1.99 ±0.23 


12.2 ±0.7 


SMG02 


22 ±3 


12.0 ±0.2 


1.41 ± 0.33 


11.9 ±0.1 


SMG04 


29 ± 10 


12.0 ±0.6 


1.80 ±0.74 


11.53 ±0.02 



hole masses for SMG06 and SMG1 1 of log(M BH /M ) = 8.7 and 
9. 1 respectively. Error s on these are taken to be equal to the intrinsic 
scatter (1.1 dex) in the lLacv et all ( 1200 ll) relationship. 

At the same time, for the non-detections (SMG10, SMG01, 
SMG02, and SMG04), we can utilise the observed correlation to 
derive upper limits on A/bh for these systems. This again requires 
we assume that no radio-quiet AGN are present, and that emission 
from large-scale jets is both negligible in these systems and in driv- 
ing the observed A/bh — Lsghz correlation, and therefore repre- 
sents a fairly rough estimate. However, keeping these assumptions 
in mind we find the 3-cr detections for SMG 10, SMG01, SMG02 
and SMG04 imply upper limits of log(M BH /M ) <: 8.1, 8.2, 7.5 
and 7.6, respectively. 

Stellar masses and their associated uncertainties (Table [5j 
were then estimated from model fits to t he observed optic al-IR 
photometry following the methodology of iDve et al.l J2008h . The 
resulting correlation with SMBH mass is pre sented in Fig.[4] alon 
with the empirical relationship determined bv lHaring & Rix feOO' 
for a sample of 30 nearby galaxies. While there is considerable un- 
certainty in the inferred SMBH masses (and also in the assumed 
Eddington ratio) owing to the large intrinsic scatter in the observed 
radio luminosity-black hole mass relation, we find that the two 
SMGs with EVN detections are statistically consistent with the lo- 
cal relation to within the systematic uncertainty in the black hole 
mass estimates. 

Of the non-detections, the predicted black hole masses are 
lower than for the EVN-detected SMG s. Only SMG10 is clearly 
consistent with fhe lHaring & Rix|j2004l) local relationship although 
the uncertainties are such that two others could also be classed as 
consistent. In only one case (SMG02) does it appear likely that the 
black hole is undermassive compared to the theoretical relation. 
However, as they are upper limits, it cannot be ruled out that all are 
underm assive, a result that wo uld be consistent with the conclu- 
sions of lAlexander et alj J2008D using X-ray phot ometry, and with 
the expectati ons of current theoretical models (Narayan an et al.l 
l20ldl2009j) i.e. that submillimetre selection preferentially iden- 
tifies starbursts, prior to the peak of SMBH activity in massive, 
gas-rich systems. 



6 CONCLUSIONS 

We have undertaken a program designed to uncover AGN radio 
cores in a sample of SMGs in the Lockman Hole, using very- 
long-baseline interferometry with the EVN. Our sensitive, high- 



resolution images have detected only two of our six targets, despite 
a strong selection bias in favour of compact emission. From their 
brightness temperatures, there can be little doubt that we are see- 
ing radio emission from active nuclei. For the other four SMGs we 
place upper limits on the radio flux density of any compact com- 
ponent associated with an AGN and conclude that star formation is 
probably the dominant source of their radio emission. 
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